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Abstract 

To investigate ferromagnetic semiconductors and insulators, such as the famous EuO, EuS, or CrBr 3 , we propose a hy¬ 
bridized Kondo-lattice model, where, in addition to the conduction electrons, localized moments (e.g., the 4/-electrons) 
are modeled as a strongly correlated band system. The quasi-empty conduction band is weakly filled due to the hybridiza¬ 
tion term. This activates the intraatomic exchange coupling between conduction and localized electrons. Temperature- 
dependent phase diagrams and quasiparticle densities of states are presented for various coupling and hybridization 
strengths. Moreover, the influence of the one-particle energy of the localized electrons Ef \s discussed. A comparison 
with mean field calculations is given at the end of this work. 
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1. Introduction 

Collective magnetism requires the existence of perma¬ 
nent magnetic moments, which arise from either local¬ 
ized or itinerant electrons. The distinction of the mag¬ 
netic moment’s type is of enormous conceptual impor¬ 
tance: While itinerant electrons cause band magnetism, 
particularly known from the classical ferromagnets Fe, Co 
and Ni, and usually examined within the Hubbard model 
Dma, local moment magnetism can be found in insu¬ 
lators and semiconductors (such as EuO, EuS and CrBra 
ISEIIS]), in the huge class of the heavy fermion systems, 
in diluted magnetic semiconductors (which are promising 
candidates for spintronic applications [71 n m HO]), and 
also in some metallic systems like Gd. Such local moment 
systems are usually described by spin models (e.g. Heisen¬ 
berg or Ising model mm) with coupling constants due to 
direct exchange or any type of superexchange. It remains 
an important aspect of modern magnetism fundamental 
research to find a unified theory which depicts the variety 
of magnetic phenomena. A small step along this path will 
be presented in this paper. 

Since collective magnetism is exclusively realized in 
solid state bodies, whose electronic structure is represented 
by energy bands and gaps that consequently determine 
the magnetic properties, it seems to be the more natural 
choice to use band models instead of pure spin models. In 
the present work, we propose a band model that describes 
semiconductors and insulators, i.e., local moment systems. 
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For our purposes, the well-known Kondo-lattice model 
[T5] . with local spins replaced by correlated /-orbitals, 
is a formidable starting point as the interaction between 
subsystems (in our case: localized and itinerant electrons) 
may lead to magnetic ordering. Note, however, that the 
interaction is of indirect nature (e.g. Ruderman-Kittel- 
Kasuya-Yosida interaction muni US], hereafter referred 
to as RKKY), thus, in search of spin-ordering effects, re¬ 
quiring at least a minimum of electrons in each subsystem. 
The problem one is now confronted with is that semicon¬ 
ductors, and insulators likewise, have quasi-empty conduc¬ 
tion bands, making indirect electronic exchange impossi¬ 
ble. In our work, we therefore extend the Kondo-lattice 
model by a hybridization term, which gives localized elec¬ 
trons the opportunity to virtually transform into conduc¬ 
tion band electrons and vice versa. This solves the problem 
of beforehand absent conduction electrons, hence allowing 
a virtual RKKY interaction, but also destabilizes the mag¬ 
netic moments in the different subsystems. As we can read 
off from our results, the last fact is of great importance for 
the understanding of the magnetism found with our model 
proposed. 

The setup of this paper is as follows: In Sec. we 
present the Hamiltonian, which characterizes the many- 
body problem that we solve by Green’s functions methods 
in Sec. Within this solution, a set of self-energies is re¬ 
quired, but so far undetermined. Approximations to com¬ 
pute the self-energies are presented in Sec. finalizing 
the theoretical framework. The results and a correspond¬ 
ing discussion are presented in Sec. A summary is given 
in Sec. El 


Preprint submitted to Elsevier 


January 6, 2015 





2. Model 


The last part of the Hamiltonian in Eq. 


For the above described systems we introduce the fol¬ 
lowing Hamiltonian H: 


Hv = Vj2[clf.^ + flc.^) , (7) 

2 ,cr 


H = Hs + Hf+ Hf^u) + Hsf + Hv ■ ( 1 ) 

The kinetic energy of the conduction band electrons is con¬ 
cerned by Hs, 


Hs — {Tij cl^Cja , ( 2 ) 

where (cia) are creation (annihilation) operators of con¬ 
duction band (s-)electrons at lattice site Ri with spin cr, 
respectively. The chemical potential is denoted by /r, and 
Tij are the usual hopping matrix elements, which are re¬ 
lated to the Bloch energies e (fc) via 

= ■ (3) 

k 

For low band occupations, a Coulomb interaction between 
s-electrons can be neglected safely. 

Assuming a non-degenerate /-level Ef, the one-partic¬ 
le energy of the localized /-electrons with creation (anni¬ 
hilation) operators (fia) is represented by 


allows for (virtual) electronic transitions from the 4/-level 
into the conduction band and vice versa, where V is the 
hybridization strength. While the s-/-exchange, Eq. 0, 
originates from the non-classical part of the Coulomb in¬ 
teraction between the conduction electrons and the local¬ 
ized /-electrons, the hybridization term, Eq. Q, is a one- 
particle scattering term which mimics, in the most sim¬ 
ple way, the hybridization between the respective bands. 
While concrete values might be obtained by standard self- 
consistent band-structure or constraint RPA calculations, 
there is a priori no reason why those couplings should not 
co-exist. 

As a consequence of the hybridization, the respective 
average occupation numbers 



a 


nfa = , 

«/ = E > 


( 8 ) 


are not necessarily constants when varying a model pa¬ 
rameter. Contrary, the total occupation number 


(4) 

2,<7 

In order to disfavor double occupancies, a large intra- 
atomic Coulomb interaction is taken into account by 

^f{U) ~ ^ ^/zo-^/i-CT j (5) 

2,(7 

where njiu = is the number operator. Since the 4/- 

wave functions’ overlap is negligible, a direct /-electron 
exchange is fairly small and hence not part of the Hamil¬ 
tonian 0. Note, moreover, that Eqs. 0 and 0 are only 
valid for systems with a total spin of 1/2, thus being an ap¬ 
proximation in modeling the 4/-levels. This simplification 
is done in order to keep mathematics on a tractable level, 
whereas the underlying physical effects are not expected 
to be affected significantly. 

To allow for magnetic ordering, interactions between 
conduction band and /-electrons must be considered. The 
s-/-exchange is commonly described by an intraatomic 
spin-spin-coupling, 

Hsf = -j'^cTi - Si , ( 6 ) 

i 

where cr^ and Si represent the conduction and /-electron 
spin operators at lattice site Ri, respectively. The coupling 
strength is given by J, where positive (negative) values of 
J stand for a preferred (anti) parallel alignment of s- and 
/-electron magnetic moments. 


n = ns+nf (9) 

is kept fixed at a constant value by proper adjustment of 
the chemical potential Since needed for later purposes, 
we also introduce the dimensionless magnetizations rugj, 

" ( 10 ) 

(7 

using Zg. = Sat - 

To fix some of the model parameters, we set the on-site 
energies of the conduction electrons to Ta = 0 eV. Hence, 
the conduction band center of gravity defines the energy 
zero. Energy units are essentially fixed by choosing W = 1 
eV for the width of the conduction band. Furthermore, we 
assume a strong Hubbard interaction Uf ^ W. 

The one-particle energy of the 4/-levels, Ef, is a pa¬ 
rameter that decisively affects the physics of the model. A 
stable local magnetic moment on the /-levels is formed in 
the limit Ef 0 and Ef + Uf ^ 0 where each /-level is 
exactly occupied by one electron. Here, however, we also 
consider a parameter regime where Ef comes close to the 
lower edge of the conduction band. This implies the pres¬ 
ence of charge fluctuations and thereby the hybridization 
term, Eq. 0 is activated which eventually generates an ef¬ 
fective /-/ magnetic exchange. Note that this implies that 
the Schrieffer-Wolff transformation HZ! does not apply to 
this parameter regime. Consequently, Hf + Hf(jj) + Ely 


2 


cannot be replaced by a local antiferromagnetic exchange 
which would trivially compete (or cooperate) with Hgf. 

It is worth mentioning that the model Q reduces to the 
conventional periodic Anderson model if the s-/-exchange, 
Eq. Q, is neglected (J = 0). On the other hand, for 
V = 0, the model essentially reduces to a fermionized vari¬ 
ant of the famous Kondo-lattice model. 


3. Theory 


In order to determine Usj and irisj self-consistently, 
we attack the many-body problem posed by the Hamilto¬ 
nian Q using Green’s function methods. With the defini¬ 
tion of the pure one-particle s- and /-Green’s functions, 


GkAE) = 




FUE) = 


E 


( 11 ) 


as well as with two mixed Green’s functions, 


( 12 ) 


all given in fc-space with N being the number of lattice 
sites, one easily obtains the following equation of motion: 

O^AE) =5 (fi + (e ik) - p) G^{E) + Vpj£'{E)+ 




ik{Ri — Rj 




(13) 


In Eq. (13), the angular brackets represent a higher 
Green’s function. To find a solution of the many-body 
problem, we use the following ansatz according to the rig¬ 
orous Dyson equation [IB] : 


[Cia, E[sf]_ )■ {E)Cma > (14) 


which is only valid within the expression for the Green’s 
function (J^Cia,Hs{\_ . The equation above de¬ 

fines the electronic self-energy E-^^^(A). Note, that Eq. 
(141 does not represent an operator identity. Inserting the 
self-energy ansatz (14) into the equation of motion (13), 
one obtains: 

Gk.{E) =Gl{E) + ^-GI{E)Y.^:Y{E)GUE)+ 


+ ^Gl{E)P^^J{E) 


(15) 


=G 




where is the Fourier-transformed self-energy, 

Gl is the correlation-free s-Green’s function. 


gUe) = 


E + ^ — e{k) 


(16) 


and where G^~° stands for the hybridization-free {V = 0) 
conduction band electrons’ Green’s function: 


G 


v=o 

ka 


(E) 


h 

E + y.-e{k)-i:^^YiE) ' 


(17) 


Since needed for later purposes, we also introduce the free 
propagator Go{E): 


= • (18) 
k 


In analogy to the system of conduction band electrons, 
the equation of motion of the localized electrons can be 
derived: 

FkAE) =^(h+{Ef- m) FUE) + VPil\E) + 


Jl ^ ^ik{Ri-Rj) ^ 



A)) J 


(19) 


As in Eq. (14), we use the following substitutions for the 
/-electrons: 


[/*a,77/((7)]_ -^ 


m 


( 20 ) 


which are only valid within the Green’s functions. Defining 
a hybridization-free /-Green’s function , 


pV=0 
^ k<7 


{E) = 




, ( 21 ) 


with {E) and y\^Y^{E) being the Fourier-transfor¬ 

med expressions of the self-energies given through Eq. 
(20), the full /-Green’s function can be written as 


F^E) = F^-AE) • 1 + 


vpI^}{e) 


( 22 ) 


The mixed Green’s functions P^} and P^^J are yet un¬ 
known. Using the substitutions (14) and (20) in the equa¬ 
tions of motion for pj^^J and EnW and using the hybri¬ 


dization-free Green’s functions (|17| and ( |21[ ), one gets: 

FkAE) , (23) 


pL\e) = 


VGl=AE) 
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and 

pUHe) = . (24) 


No higher Green’s functions appear in these results, which 
enables us to solve the system of s- and /-Green’s function 
equations by inserting Eqs. (23) and (24) into Eqs. (15) 
and (22): 


ansatz (14), the following self-energy was obtained: 


„(s),sf 

^ka 


(E) 


J 

-ZaTO/ + 

J2 Sf{Sf + 1) - + 1) 

(28) 


GkAE) = [av^) - ’ (25) 

EkAE) = . (26) 

Of course, this represents a reformulation of the full 
many-body problem only, since expressions for the self¬ 
energies or rather the hybridization-free Green’s functions 
and are still missing. We will work them out 

in the next section. 


4. Self-energies 


The last sections’ preparatory work shifts the focus 
from the original problem towards the set of (so far un¬ 
determined) self-energies, which have to be calculated ap¬ 
proximately. 

To find an analytical solution, we make an effective 
medium approach, meaning that the self-energies are 
worked out within the two respective subsystems, assum¬ 
ing E = 0. The main idea is that hybridization effects 
are already taken into account for the Green’s functions 


to some extent (see Eqs. (25) and (26)). Physically speak¬ 


ing, conduction band electrons then encounter an effective 
/-electron potential and vice versa. 

For the hybridization-free conduction band system, 
which is described by the Hamiltonian parts and (§, 
the associated Green’s function is given by Eq. ( [T7| ). Us¬ 
ing the static mean field (ME) approximation, one obtains 
the following expression for the s-electron’s self-energy: 

5.W.sf(^) ^ ^ ^ , (27) 


where stands for the z-projection of the /-system’s 
spin operator Sf. In this paper, the mean field result is 
introduced for comparison only. 

A more reliable expression for was given by 

Nolting et. al |19j . who found an approximate solution to 
the self-energy by interpolating between several exactly 
known cases, such as the zero-bandwidth limit, the fer- 
romagnetically saturated semiconductor and the second- 
order perturbation theory. The derived self-energy is trust¬ 
worthy for arbitrary temperatures and coupling strengths 
J, as well as for low band occupations rig. Likewise using 


Here, Go is given by Eq. (18) and Sf is the /-spin quan¬ 
tum number. As mentioned in Sec. we treat the /-band 
as a spin-l/2-system, i.e. Sf = 1/2. Note, however, that 
Eq. is not restricted to this case. 

The next step is finding the self-energies (E) and 
^ka'^^E) of the hybridization-free local moment system. 
Eq. ( [2T| ) shows, that knowing the /-electron’s Green’s 
function is sufficient for obtaining the sum of the wanted 
self-energies. The latter are approximated by the self¬ 
energies of the fermionized Kondo-lattice model in the 
zero-bandwidth limit. The resulting self-energy is k- 
independent. This approximation is exact for the case of 
an empty conduction band and should remain reliable for 
low conduction band filling. 

Starting once again with a mean field (ME) treatment 
of the local moment system, one obtains a (fe-inde-pendent) 
/-electron’s Green’s function which is equivalent to the so¬ 
lution of the zero-bandwidth Hubbard model [THl : 


\E + /J,-Ef+ f Zo-TTls 

^ S (29) 


E + - Ef + {zcrTTls - Uf j 

As mentioned before, this result will be primarily impor¬ 
tant for later comparison with the effective medium ap¬ 
proach. The theory of the latter, of course, also requires a 
local moment Green’s function, on which we will focus in 
the following paragraphs. 

The system turns out to be exactly solvable in the zero- 
bandwidth limit [5D]. This yields the following Green’s 
function: 

fE=^ {E) = Fj=^ (E) = ■ ( 30 ) 

Coulomb- and s-/-interaction are responsible for a split¬ 
ting into four sublevels, which are located at energy posi¬ 
tions = 

El =To - M - ^ • 5c , 

E 2 =To ~ M + w • {Sc + 1) I 

2 ( 31 ) 

E 3 =Tq — ij- + U — — ■ {Sc + 1) , 

Ei =Tq — fl + U + — • Sc , 
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each of them connected to a spectral weight ajcr, where 
j = l,...,4: 

1 


C^lcr — 


2Sc + 1 
1 


Q!2cr — ; 


(s, + l + ^ + A.^-{S, 


l)n./_cr^ , 


2Sc + 1 

*^ 3 o- “25 IjTY ~ A—ct) , 

‘^4cr = 2^ _j_ ^ (A-o- + {Sc + l)''4/-o-) • 

(32) 

In the equations above, Sc = 1/2 stands for the total 
spin of the conduction band system. The average occu¬ 
pation numbers nf-^ and conduction band magnetization 
rus were already defined in Eqs. ® and (10), whereas the 


parameter A^- is a higher correlation function: 

A, = (^S:flj,)+z,{S:n,) . (33) 

Here, S^ = S^ + iz^S^ is the ladder operator for a con¬ 
duction electron’s spin Sc, and S^ its z-component. In 
the conventional Kondo-lattice model, this two-particle 
correlation can, surprisingly, be expressed via the single 
/-electron’s Green’s function pi] : 


A„ = 


1 


JnN 


00 

^ J f.{E){E-Ef) lu, Fka{E - ii)dE , 


(34) 

where f-{E) is the Fermi function. By combining Eq. 
^0| with Eq. (34), and after tedious, but straightforward 


work, one obtains Ag-, as a functional of n/o-, risa (and thus 
ms) and A„ itself. Instead of calculating Usja within the 
hybridization-free subsystems, we chose to evaluate them 
in the full system in order to allow for a better inclusion 
of hybridization effects. We therefore use the respective 
Green’s functions (25) and ( p^ , as well as the spectral 
theorem m- 


'^S(7 - 


f-{^)^'Yl Gka{E-fi)dE 


(35) 


nfa=-^ f f-{E) Im E^{E - fi)dE . 


The fe-summation in Eq. (35) can be replaced by an inte¬ 


gration over X, weighted by the free Bloch density of states 
00 - 


1 


0o{x) = — ^5(a:-e(fc)) . 


(36) 


In self-consistent, numerical calculations, the occupa¬ 
tion numbers ngja (and hence the magnetizations, see Eq. 
( |10[ )) can now be evaluated. A non-vanishing total mag¬ 
netization, 

mtot = ms + mf , (37) 

indicates a ferromagnetic phase. 


5. Results and Discussion 

Referring to systems like EuX (X=0,S,Se,Te), a face- 
centered cubic (f.c.c.) lattice structure is considered. Cal¬ 
culations are done for an intraatomic Coulomb interaction 
among the localized electrons of G/ = 10 eV according to 
experimental data (see e.g. [H]). Recall that the conduc¬ 
tion band center of gravity is fixed at Tq = 0 eV and that 
the width of the conduction band is assumed as IT = 1 
eV. Keep also in mind, that all calculations are performed 
for a total spin of 1/2 for both conduction band (Sc) and 
local moment system (Sf). 

We divide this Section into several parts, beginning 
with metallic systems without hybridization influences 


(part 5.1). Part 5.2 addresses the case of a semiconductor 
or insulator but still T = 0. Finally, the results of the 
full theory with finite hybridization are discussed and also 


compared with mean field calculations in part 5.3 


5.1. Magnetic response of a metallic conduction band for 
V=0 

The first special case to be discussed is the metallic 
Kondo-lattice limit [H US], where the hybridization is 
switched off {V = 0). The conduction band is taken as 
half-filled. 

Some characteristic features of the half-filled Kondo- 
lattice model can be understood by looking at the quasi¬ 
particle density of states shown in Fig. For simplicity, 
we assume to have a fixed /-electron’s magnetization to/, 
which - instead of being calculated self-consistently - is 
treated here as an external parameter. Its temperature 
dependence is given via a Brillouin function with a Gurie 
temperature of 300 K. 

The spin resolved quasiparticle densities of states 
shown in Fig. split into two subbands which are roughly 
located at energies —JSf/2 and J{Sf + l)/2. This is due 
to the interpolating self-energy ansatz Eq. (28) and is 


reminiscent of the splitting found in the zero-bandwidth 


limit. A known shortcoming of the self-energy Eq. (28) is 


that even for weak coupling strengths J the gap between 
both subbands persists (see Fig. |l(a)[ ). Nevertheless, this 
does not seriously affect our weak-J calculations since the 
gap is negligibly narrow. 

By comparing the results for different to/, one finds 
the expected proportionality to the conduction band’s po¬ 
larization TOg. Note, that decreasing to/ goes along with 
an enhancement of temperature. Moreover, the spectral 
weight of the upper f-band decreases as the temperature 
is reduced, vanishing in the case of T = 0 K. This effect 
is similar to the special case of the ferromagnetically satu¬ 
rated semiconductor (for further details see Refs. [Ml 1121) ■ 
In such saturated systems (to/ = I), the /-electron cannot 
flip its spin, and thus the lower f-quasiparticle density of 
states has the shape of the free (f.c.c.)-Bloch density of 
states. Of course, this is only valid for the saturated elec¬ 
tron sort. 
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(a) J = 0.1 eV 



energy [eV] 


(b) J = 1 eV 

Figure 1: Spin-resolved quasiparticle densities of states (QDOS) of 
the conduction electrons in the Kondo-lattice model for a half-filled 
conduction band and for different coupling strengths J with fixed 
/-magnetizations rrif = 0.5 (solid black lines) and my = 1 (dotted 
red lines). The chemical potential /i is marked by the dashed vertical 
line, /i is only very slightly temperature—dependent, i.e. dependent 
on my. 


The quasiparticle density of states changes qualitati¬ 
vely in the weak coupling regime, since both subbands 
(almost) merge (see Fig. 1(a)). When increasing the cou¬ 
pling strength J, a correlation gap occurs, which separates 
the upper band from the Fermi edge (see Fig. |l(b)[ ). Then, 
exclusively the lower band determines the electronic prop¬ 
erties, which results in a rather constant magnetization 
rUg. Later we will see, that this saturation effect plays an 
important role in our model. 

These observations have been made for a given /-mag¬ 
netization and illustrate the reaction to an effective exter¬ 
nal magnetic field. Of course, also the back reaction of 
the s-magnetization on the local moment system must be 
considered in addition to get a self-consistent picture. 


5.2. Self-consistent treatment of the Kondo-lattice limit 
(V=0) 

Let us now come back to the full theory for V = 0. 
Here, ruf is no longer treated as an external parameter. 
Both TO/ and mg are rather obtained from the subsystem 
occupation numbers which are calculated self-consistently. 
To address the case of a semiconductor or insulator, the 
total occupation number is fixed at n = 1 by adjusting the 
chemical potential accordingly. Let us discuss two param¬ 
eter regimes: 

If Ef ^ /i ^ Ef + U, there are well-defined local 
magnetic moments formed in the /-electron system. This, 


however, implies that n/ = 1 and that the conduction 
band is empty. A ferromagnetic state cannot emerge in 
this case. 

Ferromagnetism is possible for V = 0 and n = 1, in 
principle, if there is a non-vanishing energy overlap of 
the lower /-electron subband with the (split) conduction 
bands. Otherwise, conduction band electrons are missing 
to mediate an effective magnetic exchange between the lo¬ 
cal moments in the /-subsystem. As the total filling is 
fixed at n = 1, the total magnetization mg + nif < 1. On 
the other hand, an overlap between s- and /-type bands 
not only implies Ug > 0 but at the same time n/ < 1. This 
means that there are charge fluctuations which destabilize 
the formation of /-magnetic moments and therefore dis¬ 
favor ferromagnetic order. Concluding, we expect that 
ferromagnetic order is unlikely for V = 0 (and n = 1), 
irrespective of the /-level position Ef. 

To address this question numerically, we have studied 
the model in an Ef-range where ferromagnetism could in 
principle be possible according to the discussion above. 
We have scanned the entire regime of coupling strength 
J. If at all, ferromagnetism was expected for strong J. 
However, only paramagnetism has been found, irrespective 
of the value of the coupling strength J. From the numerical 
point of view, this result crucially depends on the fact that, 
unlike in similar approaches (e.g. [25)1. our local moment 
system is treated fully self-consistently. 

We conclude that the self-consistent approach well de¬ 
scribes the essential spin-flip processes preventing a fer¬ 
romagnetic order of local moment and conduction band 
electrons. 

5.3. General case (V > 0) 

As the discussion above has shown, a hybridization- 
free system cannot explain the ferromagnetic nature of 
magnetic semiconductors and insulators, since there are 
no charge carriers within the conduction band (due to 
the finite gap between the latter and the local moment 
energy positions). By switching on the hybridization- 
term, hybridization-induced virtual transitions of local¬ 
ized /-electrons into the conduction band lead to a non¬ 
vanishing amount of charge carriers that participate in 
s-/-interaction processes. These are responsible for the 
alignment of the electron spins. 

The numerical studies presented below are based on an 
implementation of the full theory, as detailed in Secs. 
and|^ We again fix the total occupation number at n = 1, 
but systematically vary the energy position of the /-level 
Ef as well as the coupling strength J and the hybridization 
V. For the Hubbard-C/, we choose 1/ = 10 eV as above. 

A typical example for the resulting quasiparticle den¬ 
sities of states is shown in Fig. For two different values 
of J, referring to the weak- and the intermediate coupling 
regime, a ferromagnetic solution is easily stabilized. This 
is achieved with a small but finite hybridization V = 0.3 
eV. Note, that with Ef = —1 eV there would not be any 
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Figure 2: Quasiparticle densities of states (QDOS) in the full system 
for different coupling strengths J with fixed hybridization V = 0.3 
eV. The local moment energy lies at Ef = —1 eV. Coulomb repulsion: 
U = 10 eV. Note that /—levels located at Ef + U are not shown. 
The chemical potential fi is located at /i = —1.074 eV for J = 0.2 
eV (upper figure) and /r = —1.276 eV for J = 1 eV (lower figure). 


(b) mtot for fixed hybridization V = 0.3 eV 

Figure 3: Total magnetization mtot = m^ + mj^, represented by the 
color scale, for various s-/-interaction strengths J and hybridiza¬ 
tions V. Ej is fixed at -1 eV and the total occupation number is 
fixed at n = 1. FM is the ferro-, PM the paramagnetic phase, both 
separated by the black line, which marks the Curie temperature. 


energy overlap of s- and /-subsystem bands for V = 0. 
The fact, that the s- and /-quasiparticle densities of states 
have the same support is a mere consequence of the finite 
hybridization. 

Variation of the s-/-coupling strength J and the hy¬ 
bridization V has significant influences on the total mag¬ 
netization, and thus on the Curie temperature, as can be 
seen in our results in Fig. Summarizing, our main re¬ 
sults are: 


(i) Ferromagnetism occurs for finite V only. 

(ii) The Curie temperature increases with hybridization 
strength V. 

(iii) Around J = 0 eV the Curie temperature takes its 
highest values and the system is highly polarized in 
the low-temperature regime. 

(iv) For large J, the Curie temperature saturates. 

(v) Positive values of J show ferromagnetism, negative 
ones commonly paramagnetism. 


The limiting case V = 0 eV is a clearly paramagnetic 
situation, as discussed in Sec. 5.2 see Fig. |3(a)| For V > 0 
eV, virtual s- and /-electron transitions occur, causing a 
broadening of the quasiparticle densities of states, which is 
clear since the hybridization-term in the Hamiltonian 0 


resembles a hopping process. Moreover, the spin-resolved 
quasiparticle densities of states are repulsively shifted (as 
V acts like a perturbation), see Fig. The latter effect 
increases the spin asymmetry in the quasiparticle densities 
of states, and favors ferromagnetism. This explains point 
(ii) of the above list. 

Regarding the J-dependence, let us emphasize once 
more that conventional coupling mechanisms such as 
RKKY in the weak-J regime or double exchange in the 
strong-J regime do not apply here. The point is that the 
total occupation number n is fixed to 1. Therefore, any 
increase of the /-electron number implies a decrease of the 
conduction electron number. This is found to destabilize 
the magnetic order. 

For J = 0 eV, ferromagnetism is found to be very sta¬ 
ble with a high Curie temperature (see Fig. |3(b)[ ). Fer¬ 
romagnetic order originates in this case from the strongly 
correlated /-band where the Hubbard-C7 is much larger 
than the effective bandwidth resulting from the virtual 
hybridization processes introduced by the finite V. 

A finite J brings up two contrary effects, a splitting 
of t" and /-quasiparticle densities of states due to the 
linear-in-J (Ising-)term on the one hand and a reduc¬ 
tion of ordered spin structures owing to spin-flip processes 
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given by terms with quadratic J on the other hand. From 
our results we must conclude, that the latter dominates. 
As spin-flip processes close the gap between f- and |- 
quasiparticle densities of states, the system’s total polar¬ 
ization decreases. Nevertheless, the Curie temperature 
saturates for J —)■ oo and a finite but small magnetic mo¬ 
ment is retained (see point (iv)). 

Last but not least, the sign of J is quite important 
when looking for ferromagnetic phases. Starting off with 
the hybridization-induced shifts of the and 4,-quasi- 
particle densities of states, negative coupling strengths 
produce an Ising term that splits the densities of states 
into the opposite direction, hence reducing the spin asym¬ 
metry. As the sign of J is of importance only in the Ising 
term, the paramagnetic phase for J < 0 and the ferro¬ 
magnetic phase for J > 0 result from this first-order-in-J 
term only, just like static mean field calculations would 
predict. 

Another way of understanding our results is as follows: 
For a system with non-vanishing hybridization and a ma¬ 
jority of /^-electrons, a majority of occupied s^-states 
evolve due to hybridization. Such a parallel alignment of 
most electrons is either energetically favorable (J > 0 eV) 
or unfavorable (J < 0 eV). In the latter case, both electron 
sorts tend to order antiparallelly, which makes the occu¬ 
pation of states with opposite spin configuration more at¬ 
tractive. The less effective level-repulsion effect described 
above leads to a decrease of the spin asymmetry. 

Ferromagnetism in this system strongly depends on Ay, 
i.e. on the energy distance between the local moment and 
conduction band system. Closing (opening) the gap be¬ 
tween both subsystem’s densities of states increases (de¬ 
creases) hybridization influences. Hence, shifting Ay leads 
to an effective hybridization, which eventually changes the 
magnetization and the Curie temperature (see Fig. 4(a)). 
With increasing Ay the s- and /-subbands start to over¬ 
lap. On the one hand, this produces charge carriers in 
the conduction band which promote ferromagnetic order 
as they mediate an effective magnetic coupling between 
the /-electrons’ spins. On the other hand, increasing Ay 
simply introduces charge fluctuations in the /-levels and 
destroys the local /-moments. This leads to a rapid break¬ 
down of the ferromagnetic phase as soon as the /-level is 
pushed into the conduction band. A similar Ay depen¬ 
dence, though with different values of the Curie temper¬ 
ature, has been found in previous works using the mean 
field ansatz (see particularly Refs. [M1I2Z1)- 

As our model is not necessarily restricted to semicon¬ 
ductors and insulators (such as the Europium chalcogenides), 
we additionally investigate the physics for different occu¬ 
pation numbers n. As is seen in Fig. |4(b) a spectacular 
influence on the Curie temperature is found. Again, one 
needs to distinguish between the open gap situation (Ay 
below the conduction band) and the closed gap regime (Ay 
inside the conduction band). In the first case, decreasing 
n means to introduce more and more charge carriers in the 
conduction band system that mediate an effective coupling 



Ef [eV] 


(a) Fixed parameters V = 0.3 eV and n = 1. 



(b) Fixed coupling J = 0.2 eV and hybridization V = 0.3 eV. 


Figure 4: Curie temperature as a function of the local moment energy 
level position Ef for different parameters such as the s—/-coupling 
J (Fig. g a)|| and the occupation number n, Fig. |4(b)| 


of the /-electron moments. Hence, Tc increases, see Fig. 
|4(b)| Further decrease of n reduces the number of occupied 
/y-states, and will finally lower the Curie temperature. 

When pushing the /-level into the conduction band 
(second case), the ferromagnetic phase breaks down faster 
for small occupation numbers than for larger ones. This is 
due to the filling order, starting now with the conduction 
band (which in this situation gives the lowest states to 
be occupied). As soon as all s-electron states are filled, 
the remaining electrons will occupy the /-states given at 
energy Ay within the conduction band. The smaller n is 
chosen, the fewer /-states will be occupied, hence reducing 
the corresponding magnetization my. 

As a last result, we want to point out the differences 
between mean field theory (see also Refs. [11I17]) and our 
effective medium approach. As seen in Sec. our ansatz 
goes beyond mean field as we allow for higher-order con¬ 
tributions of the s-/-electronic correlation. Fig. shows 
that both theories give the same results in case of J = 0 
eV. As the self-energy vanishes in both cases, this is a 
trivial result (see Eqs. (27) and (28)). Strong deviations 
are found for J ^ 0 eV. All differences are due to the 
quadratic-in-J term contributions. The saturation effect 
for large J, for instance, is not found in the mean field 
results. The Curie temperature rather increases with in¬ 
creasing J without any limit. 

Even with the mean field treatment of the spin-spin 
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J[eV] 


Figure 5: Curie temperature as a function of the coupling strength 
J for mean field (dashed lines) and effective medium (solid lines) 
approach. Red graphs are calculated for constant hybridization V = 
0.3 eV, black lines for V = 0.5 eV. Ef is fixed at -1 eV. 

coupling, Eq. ([^, and even at J = 0, there is a non-trivial 
correlation effect that is induced hy Uf. The ferromagnetic 
order in the case J = 0 solely originates from this Hubbard 
interaction. In fact, due to the charge fluctuations intro¬ 
duced by V, the /-level develops a finite but small width 
such that one can effectively think of a strongly correlated 
Hubbard model with a small band width which apparently 
sustains a stable ferromagnetic phase. 

6. Summary 

We have presented a band model for describing fer¬ 
romagnetic semiconductors and insulators. Referring to 
local moment systems like EuO, EuS, Gd, etc., the model 
consists of two subsystems, namely a band of conduction 
(s-)electrons on the one, and strongly localized /-electron 
levels on the other hand. With respect to the above named 
class of substances, local moment ferromagnetism requires 
singly occupied /-levels, which can be achieved by use of 
sufficiently large on-site Coulomb interaction. Although 
the conduction band is rather sparsely occupied, or even 
empty in the extreme case of temperatures close to ab¬ 
solute zero, we expect it to be necessary for ferromag¬ 
netism in local moment systems. We therefore include 
a hybridization term which allows for (virtual) electron 
transformations between both subsystems. Last but not 
least, the two subsystems are given the chance to interact 
with each other (s-/-interaction). A hybridized Kondo- 
lattice model perfectly matches the requirements above 
and is, hence, used in our work. 

We solve the corresponding many-body problem by 
firstly treating both subsystems separately under exclu¬ 
sion of hybridization, and then using the latter to con¬ 
nect the two systems. Since each subsystem then resem¬ 
bles an effective influence on its counterpart, we call this 
an effective medium approach. In order to solve the sys¬ 
tem of conduction band electrons, a self-energy containing 
the s-/-interaction influences is required. We choose an 
(approximate) expression for the self-energy, taken from 


[25] , which results from an interpolation between a large 
amount of exactly solvable special cases, such as second 
order perturbation theory, the zero-bandwidth limit and 
the magnetic polaron. As for the local moment system, the 
strong localization of /-electron’s wave functions allows us 
to use an analytically exact zero-bandwidth solution. By 
connecting the subsystems with finite hybridization in our 
full theory, we are able to calculate occupation numbers, 
polarizations, correlations, etc., in a self-consistent way. 
Note in this context, that in particular magnetic ordering 
was investigated in a self-consistent manner. 

Our numerical calculations show that ferromagnetism 
can be found with the theory presented in Secs. [^and|^of 
this paper. In addition, magnetic ordering (and thus the 
Curie temperature) of systems with energy gap between 
/-level and conduction band is dominated by the chosen 
hybridization strength. Surprisingly, the s-/-interaction 
strength rather acts as a destabilizing parameter. This is 
most likely due to our self-consistent treatment of the /- 
electron’s magnetic moments, which therefore differs from 
conventional RKKY behavior. Contrary, fully closing the 
energy gap between /- and s-band increases electronic 
fluctuations and, hence, destroys the ferromagnetic order¬ 
ing. 

As our band model shows highly interesting physical 
behavior, the model could be used for further research on 
the system’s transport properties, thermodynamic quan¬ 
tities, or antiferromagnetic ordering properties. Another 
possible extension could bring different self-energies into 
focus, for example a fc-dependent one (see [Ml), or a self¬ 
energy which can be obtained via moment conserving de¬ 
coupling approach treatment. Also, inclusion of realistic 
DFT-results would improve the description of real sub¬ 
stances, as previous works have shown (see for example 

[Ml ISO]). 
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